Directed percolation criticality in turbulent liquid crystals 
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We experimentally investigate the critical behavior of a phase transition between two topologi- 
cally different turbulent states of electrohydrodynamic convection in nematic liquid crystals. The 
statistical properties of the observed spatiotemporal intermittency regimes are carefully determined, 
yielding a complete set of static critical exponents in full agreement with those defining the directed 
percolation class in (2+1) dimensions. This constitutes the first clear and comprehensive exper- 
imental evidence of an absorbing phase transition in this prominent non-equilibrium universality 
class. 
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Transitions into an absorbing state, i.e. a state from 
which a system can never escape, arise from simple mech- 
anisms expected to be widespread in nature. Examples 
abound in a wide variety of situations in physics and 
beyond [H, ranging from spreading processes like epi- 
demics and forest fires, to spatiotemporal chaos, catalytic 
reactions, and calcium dynamics in cells, etc. More- 
over, a host of problems such as synchronization self- 
organized criticality and wetting Q, can be mapped 
onto them. Having no equilibrium counterparts, absorb- 
ing phase transitions are central in the ongoing search for 
the relevant ingredients determining universality out of 
equilibrium. 

Over the past 25 years or so, it has been well estab- 
lished both in theory and in simulations that the vast 
majority of absorbing phase transitions share the same 
critical behavior, constituting the so-called directed per- 
colation (DP) universality class This is not surprising 
since the DP class corresponds to the simplest case of a 
single effective absorbing state in the absence of any extra 
symmetry or conservation law, as conjectured by Janssen 
and Grassberger 3j and demonstrated by hundreds of 
numerical models [l|. 

However, the situation is quite different at the experi- 
mental level. Searching for evidence of DP critical behav- 
ior, a string of experiments have been performed [7J, but 
the results remained unsatisfactory: they could not yield 
a complete set of critical exponent values in agreement 
with those defining the DP class (Table SI Q). In view 
of this state of affairs, the importance of finding at least 
one fully convincing experimental realization of DP-class 
scaling has been stressed jjjj]. The main difficulty prob- 
ably stems from the fact that one must exclude long- 
range interactions [l[ , work with macroscopic systems to 
tame quenched disorder floj| . and study them over long- 
enough scales so that the critical behavior is unambigu- 
ously measured. We have overcome these difficulties and 
report here on a firm experimental observation of DP 
criticality. 



We chose to work within electrohydrodynamic (EHD) 
convection regimes, which occur when a thin layer of ne- 
matic liquid crystals is subjected to an external voltage 
strong enough to trigger the Carr-Helfrich instability 11 1 . 
One advantage of this system is that very large aspect ra- 
tios can easily be realized, and that typical timescales are 
small (of the order of ten milliseconds). We focused on a 
transition between two topologically different turbulent 
states, called dynamic scattering modes 1 and 2 (DSM1 
and DSM2), observed successively as V, the amplitude of 
the voltage, is increased [111, [13]. The difference between 
the two states lies in their density of topological defects 
in the director field (Fig. QJi). In the DSM2 state, these 
defects, called disclinations, are created and elongated 
considerably due to shear, leading to the loss of macro- 
scopic nematic order and to a lower light transmittance. 
In DSM1, they are hardly elongated and their density 
remains very low. 

Our quasi two-dimensional cell is made of two par- 
allel glass plates separated by a polyester film spacer 
of thickness d = 12 fim. The inner surfaces are cov- 
ered with transparent electrodes of size 14 mm x 14 mm, 
coated with polyvinyl alcohol and then rubbed in or- 
der to orient the liquid crystals planarly in the x direc- 
tion. The cell is filled with A r -(4-methoxybenzylidene)- 
4-butylaniline (MBBA) doped with 0.01 wt.% of tetra- 
rt-butylammonium bromide, maintained at temperature 
25 °C with a standard deviation of 2 x 10 -4 °C, and illu- 
minated by a stabilized light source made of white LEDs. 
A CCD camera records the light transmitted through the 
plates. The observed region is a central rectangle of size 
1217 /mix 911 fim. Since there is a minimum linear size of 
DSM2 domains, namely d/\/2 we can roughly esti- 
mate the number of degrees of freedom to be 1650 x 1650 
for the convection area and 143 x 107 for the observa- 
tion area. Note that the meaningful figure is that of the 
total system size, which is at least four orders of mag- 
nitude larger than in earlier experimental studies Q ■ In 
the following, we vary V and fix the frequency at 250 Hz, 
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FIG. 1: (Color online) Spatiotemporal intermittency between DSM1 and DSM2. (a) Sketch of a DSM2 with its many entangled 
disclinations, i.e. loops of singularities in orientations of liquid crystals, (b) Snapshot taken at 35.153 V. Active (DSM2) regions 
appear darker than the absorbing DSM1 background. See also Movie SI [g]. (c) Binarized image of (b). See also Movie S2 [jj]. 
(d) Sketch of the dynamics: DSM2 domains (gray) stochastically contaminate [c] neighboring DSM1 regions (white) and/or 
relax [r] into the DSM1 state, but do not nucleate spontaneously within DSM1 regions (DSM1 is absorbing), (e) Spatiotemporal 
binarized diagrams showing DSM2 regions for three voltages near the critical point, namely 34.858, 34.876, and 34.900 V. The 
diagrams are shown in the range of 1206 fj,m x 899 fjxa (the whole observation area) in space and 6.6 sec in time. 



roughly one third of the cutoff frequency 820 ± 70 Hz. 

Closely above the threshold V c marking the appear- 
ance of DSM2, a regime of spatiotemporal intermittency 
(STI) is observed, with DSM2 patches moving around in 
a DSM1 background (Fig. [lb and Movie SI [§]). This 
suggests an absorbing phase transition [3] with DSM1 
playing the role of the absorbing state (Fig. [Hi) . Prior to 
any analysis, we must distinguish DSM2 domains from 
DSM1. This binary reduction can be easily performed by 
our eyes, so we automated it using the facts that DSM2 
regions look darker, have l ong er time correlation, and 
have minimum area of <i 2 /2 [l2| (see Ref. Q for details). 
A typical result is shown in Fig.[TJ; and Movie S2 Fig- 
ure [1^ shows spatiotemporal diagrams obtained by that 
means in the steady STI state. They illustrate that the 
key ingredient is present: active DSM2 patches evolve 
in space-time essentially by contamination of the DSM1 
absorbing state and recession from the DSM2 state (as 
sketched in Fig. |TJi) . In this context, the order param- 
eter p is just the ratio of the surface occupied by active 
(DSM2) regions to the whole area. 

We first observe the steady-state STI regime under 
constant voltages V, in the range of 34.858 V < V < 
39.998 V. The voltage for the onset of steady roll convec- 
tion is V* — 8.95 V. The time-average of p in the steady 
state, p, is taken over the period 1000 sec < T < 8000 sec, 
which is longer than 1000 correlation times. It shows a 
continuous and algebraic decay to zero with decreasing 
V, a clear signature of criticality (Fig. [2]). The critical 
voltage V c and the critical exponent (3 are determined 
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FIG. 2: (Color online) Variation of the average DSM2 fraction 
p with V in the steady state. Inset: same data in logarithmic 
scales, with vertical errorbars indicating the standard devia- 
tion of the time series p(t). Fluctuations of V are negligible 
except for the first data point, where the standard deviation is 
shown by a horizontal errorbar. Blue lines are fitting curves. 



by fitting these data to the usual form p ~ (V 2 — V 2 )^ 3 . 
(Here, deviations from criticality are measured in terms 
of V 2 instead of V by convention, since the dielectric 
torque that drives the convection is proportional to V 2 
ll|.) This yields the following estimates 



V c = 34.856(4) V, = 0.59(4), 



(1) 



where the numbers in parentheses indicate the uncer- 
tainty in the last digits, which correspond here to a 95% 
confidence interval in the sense of Student's t. Our esti- 
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FIG. 3: (Color online) Histograms of inactive (DSM1) length i x , i a and duration r in the steady state near criticality. Dashed 
lines indicate the estimated algebraic decay at threshold. 



mate (3 = 0.59(4) is in good agreement with the (2+1)- 
dimensional DP value /3 DP = 0.583(3) Q. 

We then measure N(l) and N(r), the distributions of 
the sizes / and durations r of the inactive (DSM1) re- 
gions. For instance the distribution N(l x ) in the x di- 
rection is obtained by detecting all inactive segments in 
the x direction for all y and t. We find that they de- 
cay algebraically at criticality up to a finite-size expo- 
nential cutoff (Figs.Ek-c). The power-law decay is fitted 
as N(l) ~ l-^ and N(t) ~ with 

p x = 1.10(22), p y = 1.23(4), M|| - 1-61(6), (2) 

where p x and p y indicate the exponent p± measured in 
the x and y direction, respectively. They are again in 
good agreement with the DP values f/[ F = 1.204(2) and 
//P p = 1.5495(10) [14]. These exponents are related to 
the more conventional correlation length and time expo- 
nents v± and v\\ via p± = 2 — (i/v±_ and p\\ = 2 — f3/v\\, 
which lead us to estimate 



v x = 0.66(17), v y = 0.77(7), = 1.51(25). 
They should be compared to__the DP values 



(3) 



DP 



0.733(3) and vf F = 1.295(6) [14|. The above estimates 



for f± arc in reasonably good agreement with the DP 
value, while the agreement on v\\ is less satisfactory. Our 
second set of experiments provides another, independent 
estimate of v\\. 

Setting the voltage to 60 V, i.e. much higher than V c , 
we wait until the system is totally invaded by DSM2 do- 
mains. We then suddenly decrease V to a value in the 
range of 34.86 V < V < 35.16 V, i.e. near V c , and observe 
the time decay of activity for 900 seconds. We repeat this 
10 times for each V value and average the results over this 
ensemble. In such "critical-quench experiments," correla- 
tion lengths and times grow in time, and, as long as they 
remain much smaller than the system size, the scaling 
estimates are free from finite-size effects. As expected, 
p(t) decays exponentially with a certain correlation time 
for V < 35.02V, algebraically over the whole observation 
time for V — 35.04V, and converges to a finite value for 



V > 35.06V (Fig. |4^). A simple scaling Ansatz implies 
the following functional form for p(t) in this case: 

p(<) ~r a /( £ ^ii), 



(4) 



where e — (V 2 — V 2 )/V 2 is the deviation from criticality 
and /(C) is a universal scaling function. From the slopes 
of the algebraic regimes for the three V values closest to 
the threshold, we estimate 



V c = 35.04(2) V, a = 0.50 



+(8) 

-(5) - 



(5) 



Note that V c is slightly higher than in the steady-state 
experiments. In fact, the roll convection onset V* = 
8.96 V is also higher. We believe this is because, during 
the few days which separated the two sets of experiments, 
our sample has aged, a well-known property of MBBA. 
On the other hand, no measurable shift was detected 
during a given set of experiments. Our a value is again 
in g ood agreement with the DP value a DP = 0.4505(10) 
14]. Finally, our direct estimates of (3 and a in Eqs. ([1]) 
and ((5j), and the scaling relation (j4j yield 



1.18 



+ (14) 
"(21)' 



(6) 



which is in good agreement with the DP value i/P = 
1.295(6) |l3 |. Note that our two estimates of z/y lie on 
opposite sides of the DP value, with their confidence in- 
tervals overlapping this reference number. Furthermore, 
Eq. (j4|) implies that the time series p(t) for different volt- 
ages collapse on the universal function /(£) when p(t)t a 
is plotted as a function of t\e\ v w . Our data do collapse 
reasonably well on the universal function of DP (Fig. 2b) • 
Continuous absorbing phase transitions are character- 
ized by three independent (static) exponents [l[ . So far, 
we have found those three independent algebraic scaling 
laws (typically over two decades). The critical exponent 
values all agree within a few percent with those of the 
DP class in (2+1) dimensions. In addition, data collapse 
is satisfactorily achieved onto the universal scaling. This 
constitutes the first complete and unambiguous experi- 
mental realization of a DP-class absorbing phase transi- 
tion 15]. Looking back at previous attempts to exhibit a 
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FIG. 4: (Color online) Critical-quench experiments, (a) De- 
cay of p(t) after the quench, for V = 34.86 V, 34.88 V, 
35.16 V from the bottom left to the top right. The data for 
V — 35.04 V (showing the longest scaling) are indicated by a 
thicker line, (b) Scaling plot of data in (a), with V c , a and v\\ 
values estimated from the experiment [Eqs. ([5]) and (|6]l]. The 
dashed curve indicates the DP universal scaling function /(£) 
obtained numerically from the process sketched in Fig. [T}1 
(so-called contact process). A collapse of similar quality is 
obtained when using DP-class exponent values. 



DP phase transition experimentally, one may wonder why 
we obtained such clear DP scaling laws compared with 
earlier experiments (see Table SI |8|). A major reason, al- 
ready mentioned, is that EHD convection offers effective 
sizes orders of magnitude larger than in previously stud- 
ied systems such as Rayleigh-Benard convection. Simi- 
larly, typical timescales are much shorter. But beyond 
those obvious facts, we believe that our unusual choice 
of a transition between two fluctuating states has also 
helped: in many past experiments, the absorbing state 
was essentially fluctuation- free (laminar) 7], and thus 
may have caused long-range effects (mean flows, rigidity 
of the laminar pattern), which probably further reduced 
the effective size of the systems and could even break DP 
universality In our system, such long-range interac- 
tions are likely to be killed by local turbulent fluctuations 
of the absorbing state. In short, our experimental setup 
easily reveals DP scaling because it is able to tame both 
quenched disorder (being macroscopic) and long-range 
interactions (being turbulent). 
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